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C^i Abstract 

Oh! 

r~| ! A new perturbative scheme is proposed for the evaluation of the free energy 

density of field theories at finite temperature. The screened loop expansion 
takes into account exactly the phenomenon of screening in thermal propaga- 
tors. The approach is tested in the A^-component scalar field theory at 2-loop 

C^ '_ level and also at 3-loop in the large A^ limit. The perturbative series gener- 

ated by the screened loop expansion shows much better numerical convergence 
than previous expansions generated in powers of the quartic coupling. 



PACS: ll.lO.Wx, 12.38.Mh 

Keywords: QCD, scalar field theory, free energy density, Debye screening, loop 

expansion, gap equation, high-T expansion 



1 Introduction 



The weak coupling expansion of the QCD free energy density does show very bad 
convergence properties. The coefficients of the usual perturbative series are of alter- 
nating sign and of increasing magnitude ^-^ Only in the TeV temperature range 
can one find a satisfactory numerical convergence rate 0]. This is in great contrast 
to numerical calculations of this quantity (or of its derivative, the internal energy) in 
Monte Carlo simulations. These calculations indicate that deviations from the high 
temperature ideal gas limit are within 15% already for temperatures about twice 
the critical temperature p. 

On the other hand it has been observed already quite some time ago that the 
formula for the free energy density of a massive ideal gas does give a quite satisfactory 
description of the numerical simulations 0. Such a description leads to results 
about 10% below the Monte Carlo data. Although such an approach clearly has to 
be refined in order to take care of the correct number of degrees of freedom and the 
modifications due to interactions it indicates that a representation of a perturbative 
expansion in terms of massive degrees of freedom might be a good starting point. 

This observation leads us to conjecture that the apparent poor convergence of 
the perturbatively calculated free energy density might be improved if the series is 
reorganised into a loop expansion evaluated with screened propagators; i. e. instead 
of expanding around the massless ideal gas limit, we intend to perform the loop 
expansion starting from a massive ideal gas. 

There are several screening scales in QCD, which makes the practical application 
of this conjecture more complicated, than it looks like at first sight. In the scalar 
theory, on the other hand, the only relevant mass is the Debye screening mass. For 
this reason we will discuss here only the calculation of the free energy density of the 
N component scalar $^-theory. The problem of poor numerical convergence shows 
up equally well already here [|1^, [11], |12[. One thus can test to what extent the above 
mentioned reorganisation of the perturbation series improves the convergence. 

We shall calculate and analyze first the two-loop correction to the massive ideal 
gas formula. With respect to their screening mass dependence our formulae are 
exact. The ultra-violet finiteness of our calculation is ensured without the need to 
fix the screening mass inherently. As an external information we shall use for the 



screening mass the root of the corresponding 1-loop gap equation [0, |1^, |I6| . For 
large N it also is possible to estimate the 3-loop correction of the free energy density. 
Finally, we shall compare the results obtained with this screened perturbative expan- 
sion with those obtained from the conventional resummation scheme [^, 0, |T1|, [T^ 



as well as from an effective theory approach |13 



2 Loop Expansion with Debye Screened 
Propagators 



We consider an 0{N) symmetric scalar field theory with the following Lagrangian: 

L = Lq -\- Lint, 
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Following ref. |I0| , |TT| we have introduced and substracted a mass term with m ~ 
0{g) in order to reorganize the perturbative expansion. The coupling constant 
renormalisation factor is given by |T^ 



l + ^^ + 0(.^) 
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and we have left out the field renormalisation factor Zi, which is unity up to 0{g'^). 

All Feynman diagrams contributing to the free-energy density up to 3-loops can 
be found in p[- The contribution of these diagrams are: 
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Here we have used the same labeling of the diagrams as in [Q. Furthermore, we use 
the abbreviation A^ = (A^ + 2)/3 and have introduced the following notation for the 
relevant integrals: 



J{m) 



f log(P^ + m^ 
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I'M = Ltt;:^— 7^2, (6) 
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Here we are using dimensional regularisation at finite temperature. The symbol Jp 
therefore is understood to represent T J2'^=-oo I it ~\i-'i^ ■ We note that all relevant 
sum-integrals appearing in Eqs.(3-7) are UV-singular and do depend on the renor- 
malisation scale /i. We shall argue below that a consistent expansion is obtained by 
simply dropping these terms at each order in the loop expansion. 

At two loop level only the diagrams (a-c) contribute. As will be discussed below we 
have, at present, difficulties to extend our analysis at finite N to 3-loops because we 
cannot handle the finite terms in the "basketball" diagram (h) for nonzero values of 
the screening mass. We have, however, also listed the 3-loop diagrams here, since the 
contribution from diagram (h) is suppressed in the large N limit. We are therefore 
also going to analyze the behaviour of our screened perturbative expansion in this 
limit by including the 3-loop contributions (d-g). 

We have explored two possibilities to evaluate the 1-loop integrals Jim), I{m) 
and I'{m). The first consists of performing the Matsubara sum using the contour 
trick, subtract the T = piece, and then evaluate the remaining three-dimensional 
integral numerically^. For J{m) we obtain 



> / -^^Jf^^^ + 2,^^ /^!4^1og(l-e-^^). (8) 

The other two integrals, I{m) and I'{m), can be obtained from this by simple 
differentiation with respect to m^. The other possibility to evaluate these sum- 
integrals is to use directly the high-T expansion of the integrals J{m), I{m) and 
I'im) which can be found in [ITSi. 



The leading order 1-loop result is given by diagram (a), i.e. F'^'^\ It represents 
the contribution of a massive ideal gas. However, on its own this contribution is not 
complete. It contains UV-divergent and scale dependent contributions |T^. The 



divergent piece is proportional to m^ /e. Since the screening mass is linearly depend- 
ing on T this divergent term cannot be canceled by the subtraction of appropriate 
zero temperature contributions. Higher loop contributions should take care of the 
cancellation of temperature dependent singularities. As the screening mass is 0{gT) 
we expect that two and three loop contributions are needed to cancel the leading 
poles in e. In fact, it is easy to check, that this is the case. Collecting all poles in e 



"Before performing the numerical evaluation of Eq. (8) we have to expand the integration 
measure in e. 



with residua proportional to m^, g^T'^im?, g^T^ coming from diagrams up to 3-loop, 
we find that the singular terms cancel irrespective of the actual value of m^. 

This observation forms the basis for our conjecture: The ultra-violet, temperature 
dependent singularities appearing in a given order of the screened loop expansion will 
be canceled by appropriate contributions from higher order loops. At any given order 
in the expansion these singular contributions thus can be dropped. 

At two loop order there appear also singularities with poles in e^ as well as 
poles with higher powers of m^. In order to check explicitly the cancellation of 
these singularities we would need to know the complete 4-loop contribution to the 
free energy density. This is not known. We thus invoke the above conjecture to 
eliminate these singular terms from the expansion. 

Another question which has to be investigated is the dependence on the renor- 
malisation scale. In the conventional perturbative approach the scale dependence 
appears first at 3-loop level. The scale dependent contributions appearing in the 1- 
loop and 2-loop results (see the explicit expressions of J{m) and /(m) in ||18|) usually 



are canceled by corresponding zero temperature contributions. In the screened per- 
turbative approach this again cannot be achieved because the scale dependent terms 
are proportional to the temperature, i.e. the screening mass. However, these addi- 
tional scale dependent terms, proportional to the screening mass, are canceled by 
higher loop contributions in analogy to the singular terms discussed above. Further- 
more, the screening mass independent scale dependence of the 3-loop contribution 
is exactly canceled when one takes into account the scale dependence of the cou- 
pling constant in the 2-loop diagrams. Up to C(5'^) the free energy is thus scale 



independent |13]. We therefore omit terms proportional to ln(yU /T ) in the lower 



loop expressions, tooQ. 

Finally, we have to specify the screening mass. Since m"^ has remained undeter- 
mined we have specified it as the exact root of the gap equation relevant for the 
chosen order of the loop expansion, i.e. we use the 1-loop gap equation[| without 
expanding I{m) in its argument: 

Below we shall use roots of Eq. (^ for the analysis of the free energy density. The 
screening mass obtained this way agrees well with the 1-loop value for g<l while 
it is only half as large as the 1-loop value for g ~ 10. We also note that as a 
consequence of using the root of the gap equation for m^, no term proportional to 



^/2 is the MS scale, /i^ — ^j^fl"^- 

■^To make the 1-loop gap equation self-consistent we have to substract singular and scale de- 
pendent terms also here. These terms would be generated by higher loop contributions (see 
discussion in the next section ) which are not contained in the 1-loop gap equation {m. 



g^m will appear in the free energy density. Such a term is seen to exist in F^''^ when 
we use the high-T series for I{m). However, this contribution is canceled by part 
of F'^^\ when exploiting (P). The entire cubic dependence on g^^ which partly is 
responsible for the poor convergence of the conventional perturbative expansion, is 
now hidden in a term proportional to th? which contributes to J{m). 



3 Numerical Results 

We have evaluated the free energy density using Eq. (|]) without making any as- 
sumption about the magnitude of m, i.e. we have calculated F "exactly" in m. 
For finite A^ we can do this so far only at 2-loop level. At 3-loop one would have 
to calculate the finite part of haiiifn) for arbitrary values of m. Since we are not 
assuming that m is small we did use the integral representations for J{m), I{m) 
and I'{m) in our numerical calculations. However, we also checked the alternative 
approach based on the high-T expansion and found that this series does converge 
rapidly. For the integrals listed above it is sufficient to take the first six terms of 
high-T expansion for g G [0, 12]. 

In Fig 1. we show the 1-loop and 2-loop results for the free energy density as 
function of the self-coupling g and for A^ = 1 . The screening mass has been computed 
from Eq. (H). Also shown there are the results obtained with the conventional 
perturbative expansion up to 0{g'^) and 0{g^) in the scalar self-coupling. The 
Figure clearly shows that the difference between 1-loop and 2-loop results of the 
screened perturbative expansion is small and essentially stays constant over the 
entire range of couplings. In the conventional perturbative approach on the other 
hand different orders in the expansion differ largely from each other. 

For large A^ the contribution of the basketball diagram is suppressed by 1/A^ 
relative to other diagrams. In this limit we thus can consider also the 3-loop result. 
In the large- A^ limit the 1-loop gap equation (^) becomes exact. Its derivation as 
a Schwinger-Dyson equation has been given in |T^. The setting sun diagram with 
exact vertex becomes negligible in the large N limit, while the divergent and scale 
dependent part of the tadpole contribution is canceled by the renormalisation of the 
scalar self-coupling. 

For our numerical investigation of the large-N limit we have chosen A^ = 10. 
The 1-loop, 2-loop and 3-loop results for the free energy density obtained from the 
screened perturbative expansion are shown in Fig. 2. The screening mass has been 
determined from Eq. (^. For comparison we also show the conventional 0{g^) result 
of 0. Lower orders in this approach do show alternating behaviour similar to that 
shown in Fig. 1 for the case A^ = 1. We note that like in the A^ = 1 case the 
higher loop corrections lead to smaller deviations from the massless ideal gas than 
the leading 1-loop result. We also find that for g < 12 the difference between the 
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3-loop result and the massless ideal gas limit is only 10%. 



4 Conclusions 

We have formulated a perturbative scheme for the iV-component scalar field theory 
which includes a temperature dependent screening mass in the scalar propagator. 
We conjecture that a systematic loop expansion is obtained when singular terms 
and scale dependent terms proportional to the screening mass are dropped order by 
order in the perturbative expansion. This expansion coincides to all orders with the 
conventional approach, if the screening mass is treated perturbatively. Within the 
screened perturbative approach we can handle the effect of a screening mass exactly. 

The analysis of the A^ = 1 scalar field theory and its large- A^ limit to 2-loop 
and 3-loop, respectively, shows that within this approach the convergence of the 
perturbative series is drastically improved over the conventional approach. It is 
also evident that taking into account higher orders in the massive loop expansion 
systematically brings the free energy density closer to the massless ideal gas limit. It 
certainly would be interesting to test these results against non-perturbative lattice 
calculations. 

We hope that the nice properties seen here for the scalar theory remain valid 
also in the case of QCD where more mass scales will become relevant. 
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Figure 1: Shown is the free- energy density of the N=l scalar field theory as a function 
of the scalar self- coupling g in the units of the free energy density of a massless ideal 
gas (^^^-)- The curves represent 2-loop (a) and 1-loop (h) results of the screened 
loop expansion as well as the 0{g'^) (c) and 0{g^) (d) results of the conventional 
perturbative expansion. 
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Figure 2: The free energy density of the scalar field theory in the large-N limit as a 
function of the scalar self-coupling g in units of the free energy density of a massless 
ideal gas. Shown are results for N = 10 (see text). The curves represent 3-loop (a) 
2-loop (h) and 1-loop (c) results of the screened loop expansion as well as the 0{g^) 
(d) results of the conventional perturbative expansion. 
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